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a. 

Abstract 

We derive non-linear differential equations for correlation functions of U(l) twist fields in 
the two-dimensional massive Dirac theory. Primary U(l) twist fields correspond to exponen- 
tial fields in the sine-Gordon model at the free-fermion point, and it is well-known that their 
vacuum two-point functions are determined by intcgrable differential equations. We extend 
part of this result to more general quantum states (pure or mixed) and to certain descen- 
ds j. dents, showing that some two-point functions are determined by the sinh-Gordon differential 
1} ■ equations whenever there is translation and parity invariance, and the density matrix is the 
exponential of a bilinear expression in fermions. We use methods involving Ward identities 
associated to the copy-rotation symmetry in a model with two independent, anti-commuting 
copies. Such methods were used in the context of the thermally perturbed Ising quantum 
field theory model. We show that they are applicable to the Dirac theory as well, and we 
suggest that they are likely to have a much wider applicability to free fermion models in 
general. Finally, we note that our form-factor study of descendents twist fields combined 
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, with a CFT analysis provides a new way of evaluating vacuum expectation values of primary 
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U(l) twist fields: by deriving and solving a recursion relation. 
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1 Introduction 

The correlation functions of local fields in a quantum field theory (QFT) yield all the physical 
information contained within a model. In most models these functions are only accessible 
perturbatively and the expressions obtained only valid for large energies. It is well known, 
however, that there exists a set of two-dimensional intcgrablc quantum field theories, having 
free-fermion representations, where non-trivial vacuum correlation functions can be expressed 
as solutions to intcgrablc differential equations, the first result of this nature being the two- 
point spin-spin correlation function in the thermally perturbed lattice Ising model [24] . Since 
then it has been found that this result, as well as similar results in the QFT model of Dirac 
fermions, can be retrieved through alternative methods: holonomic quantum fields (Dirac 
model) [20], Fredholm determinants (Ising and Dirac models) [1, 12, 2], determinants of 
Dirac operators (Dirac model in flat and curved space and with magnetic field) [17, 18, 13], 
and doubling of the model (Ising spin chain, and Ising QFT model at zero and non-zero 
temperature and in curved space) [19, 9, 6]. The existence of such differential equations is 
extremely useful, as correlation functions can be evaluated with very high accuracy once 
the initial conditions have been fixed using conformal perturbation theory and form-factor 
analysis [3, 6, 7, 13, 14]. Note that, at zero temperature, this provides solutions to non-trivial 
connection problems in the theory of Painleve equations. 

Thanks in great part to the holonomic quantum field description, it is understood that 
such differential equations should exist in general for correlation functions of twist fields in 
free fermion models. Twist fields are local fields which have, with respect to the fermion 
fields, non-trivial monodromy associated to a symmetry transformation. The concept of 
twist fields was first introduced in [21], as the Z2 monodromy field of the Majorana fermion 
corresponding to the spin field of the Ising model. Similarly, it is known that primary 
U(l) twist fields in the Dirac model reproduce correlation functions of exponential fields in 
the sine-Gordon model at a particular value of the coupling (the free-fermion point). It is 
correlation functions of such I2 and U(l) primary twist fields that have been studied and 
that are known to lead to differential equations in QFT models. 

In this paper, we examine the model of free massive Dirac fermions and study two-point 
correlation functions of U (1) twist fields. The method used to study the correlation functions 
is similar to that of [9] where differential equations for Ising spin-spin correlation functions 
are found by looking at conserved charges in a model containing two non-interacting copies 
of the Ising field theory. The associated Ward identities naturally give a bilinear form of 
the differential equations. However, in [9], many particularities of the Ising model were used 
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to obtain such results. The aim of this paper is four-fold: to clarify the method of [9] and 
show how it works beyond the Ising model and beyond the finite-temperature density matrix 
(suggesting that it should work more generally for any free-fermion model and quite general 
density matrices expressed as exponentials of bilinears); to establish more precisely in the 
Dirac theory differential equations for twist-field correlation functions in general quantum 
states, pure or mixed (under the condition that the density matrix is the exponential of a 
bilinear expression in the fermions); to establish equations determining two-point functions 
involving certain descendent fields; and to derive a recursion relation for vacuum expectation 
values (VEVs) of primary twist fields which can be used to evaluate them. 

The method of [9] was widely thought to be only applicable to the Ising model be- 
cause of the particular form of the OPEs: only the primary bosonic twist field, its unique 
fcrmionic descendent, and their derivatives are involved. In the Dirac theory, there are many 
more fermionic descendents, and there are bosonic descendents that are not derivatives of 
primaries. Yet, the method still works. We believe that the method should work quite 
generally because it is intimately related to standard techniques used in the context of clas- 
sical integrable models, where bilinear differential equations are naturally obtained. The 
method of [9] was also slightly restricted in other ways: it only dealt with the vacuum or the 
finite-temperature states, and it used explicitly translation invariance in order to construct 
particular conserved charges and derive differential equations from the associated Ward iden- 
tities. We show, in the context of the Dirac theory, that much more general quantum states 
can be used, and we show how to obtain equations even without translation invariance, only 
using a single conserved charge (something which was initially noticed, in the Ising model, 
in the thesis of one of the authors [4]). We note, however, that only when the quantum state 
is translation and parity invariant, these equations simplify to something that is tractable. 
Interestingly, this new way clearly separates the origin of the mass dependence on the Ward 
identities from pure CFT calculations. 

In order to give a full understanding of the way the equations are derived, we develop the 
theory of Dirac U(l) twist fields to some extent: we construct spinless, chargeless twist fields 
a a for all non-integer a € R, characterised by their U(l) monodromy e 2ma and dimension a 2 , 
as well as their fermionic descendents defined via the operator product expansion (OPE), 
c a +i,a ~ ^R a a and a a -i ta ~ \& R(j a (where is the right-moving component of the 
Dirac spinor). The construction is done both using explicit form factors, and using a CFT 
description that takes into account the constraints provided by the massive perturbation. 
We evaluate the action of the double-model conserved charge using these two pictures. 

After calculating the form factors and examining the OPEs, a recursion relation for the 
VEVs (vaclCTajvac) = m a c a is found: 



This relation is consistent with known results [15] but can also be viewed as a method for 
evaluating c a , given the normalisation cq = 1 and the symmetry c_„ = c a 

The main results are as follows. Consider correlation functions in a quantum state whose 
density matrix is an exponential of a bilinear form in fermions, and which is translation and 
parity invariant. Writing the following correlation function as 



r(a + l) 

T(-a) 



(1) 



(a a (x 7 y)a a {0,0)} 




(2) 
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we show that it can be parametrised by a single function if) where 

2 

— 777 

9dt/j = — sinh(2?A) (3) 

9 

— 777 

99S Q = — (1 - cosh(2V0). (4) 

This result was first found implicitly in [20] and then explicitly, using the method of Frcdholm 
determinants, in [2]. We further find that the correlation function involving the descendent 
field (Tq+i along with a a is obtained from 

(vacMz,yK+i(0,0)|vac) = c a c a+1 m 2a2+2a+1 e s '^ x ^ (5) 

via 

Note that this is different from the equation found for the correlation functions 
(vac|cr Q (a;, y)ap(0, 0)|vac) in [2], but there is no contradiction, as there there was the restric- 
tion < a, p < 1. Finally, we find that the correlation function of fcrmionic descendents, 

ta+i,a(z, !/K,a+i(0, 0)) =e- i ™c a c a+1 m 2a2+2a+1 n(x,y), (7) 
is algebraically expressed in terms of the previous correlation functions: 



ft = v/ e 2S Q _ e 2SJ,_ (g) 

The paper is set out as follows: in Section 2 notation and conventions are established 
and the definitions of the primary twist fields are made explicit. In Section 3 the descendent 
twist fields are introduced through the OPEs of fcrmion fields and primary twist fields. The 
one- and then multi-particle form factors are calculated and their consistency is checked. A 
recursion relation for the constants c a is found and discussed. In Section 4 the double Dirac 
model is introduced. Conserved charges are discussed and the method by which the action of 
these charges is verified is described. In Section 5 the appropriate Ward identities are written 
down and simplified to give the main result. Section 6 gives our conclusions and possible 
areas for future work. Finally the appendices contain a discussion of the CFT considerations, 
and of the Ward identities of Section 5 in a model without translation, rotation and parity 
symmetry. 



2 The model and its primary fields 

Throughout this paper, x G K is a space coordinate and y G K is a Euclidean time coordinate 
(y = it, where t is real time). We will put space and Euclidean time into complex coordinates 
z := — ^(x + iy), z := ^(x — iy), with derivatives d := d, = id x + d y and d := d z = —id x + d y . 
Also, when only one coordinate is specified for a quantum field, it is understood as the space 
coordinate, the Euclidean time coordinate being set to 0. 



2.1 Dirac fermions 

In two dimensions and in the quantisation on the line, the free Dirac Fermi field of mass m 
is an operator solution ^L(x,y) to the equations of motion 

R L (g) 
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and to the equal-time anti-commutation relations 

{^ R (x 1 ),¥ R (x 2 )}=47r5(x 1 -x 2 ), {^ L {x 1 ),^l{x 2 )}=inS(x 1 -x 2 ) (10) 

(other combinations of fields anti-commutating) , subject to the condition that there exists a 
state |vac), the vacuum state, such that 

lim ^ fljL (x,y)|vac) = lim t> R L (x, y)|vac) = 0. (11) 

y— > — oo y—^—oo ' 

The normalisation in (10) is chosen in such a way that the fields have the so-called CFT 
normalisation in terms of the z variables: as \z\ — z 2 \ — > 0, 



(vac IT 



vac|T *l(xi,tfi)*i(a:3,»3) 



zi - z 2 

1 

21 - 22 ' 



Here, the symbol T is the time-ordering symbol: in its bracket, operators at later Euclidean 
times are placed to the left of operators at earlier Euclidean times (getting a sign, in general, 
if the operators exchanged have fermionic statistics). 

As usual, the solution can be given in terms of creation and annihilation operators (Fourier 
modes of the fields), whose time evolution is simple: 

V R (x,y) = V^Jd9e s/2 (p\.(0)e vE »^ txp9 - iD-(0)e- yEa+lxpB ^ 

9 L (x,y) = V^ijd9e- e/2 (iDl(0)e yE »- ixp ° - D_(e)e~ yEs+ixp ^ (12) 

where 

Eg = mcosh#, pg^msvah.9. (13) 

It is worth noting that since y is a Euclidean time, the Hermitian conjugate operators are 
obtained from ^ Ri l by direct Hermitian conjugation of the expressions above, followed by a 
change of sign of y: 

#k(x,y) = Vmjd9e / 2 (iDl(6)e yEe - txpe + D + {d)e- yEe+lxpe ^j 

*<(*,„) - ,/*/*.-•/■ -W +m e-^ «) (14) 

From (10), the creation and annihilation operators satisfy 

{Dl(6 1 ),D + (6 2 )} = 6(6 1 -6 2 ), {Dl(6 1 ),D_(6 2 )} = 6(6 1 -6 2 ) (15) 
(with all other combinations anti-commuting), and from (11), the vacuum condition 

D± (0)| vac) = 0. 

The Hilbert space is simply the Fock space over (15), with a basis of multi-particle states 
denoted by 

\6 1 ---e n ) ei ... en :=Dt 1 (e 1 )---Dl n (0 n )\vac) for 0i > ■ • • > n , a G {+, -}. (16) 
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These states have total energies • Eg j and total momenta J2j PSj ■ The dual vectors will 
be denoted by 61 ' " e " (#i • • ■ 6 n \ := \6\ ■ ■ ■ # n )ei— e n , an d clearly we have 



n^m-^)- (17) 



Writing as in (16) vectors with different ordcrings of the creation operators, the identity 
operator can be decomposed into 

°° -i />oo />oo 

1= E1vTE/ ^■■•/ ^l^-'-Me,...^ - ei (^"--^|. (18) 

The solution (12), (14), defining the mode operators from the Fermion operators, is 
chosen in such a way that it satisfies the equations of motion, and that the resulting mode 
operators obey the canonical anti-commutation relations (15). However, the choice of such 
mode operators is not unique: we may re-define the mode operators via D±(6) H> u±D±(8) 
for some pure phases u±, \u± \ — 1, and still have a solution to the equations of motion and 
canonical anti-commutation relations. This amounts to a phase-redefinition of the states in 
the Hilbcrt space. The mode operators are partially fixed by requiring further the crossing 
symmetry condition 

(vac|* flii (O)|0 + in) t = - £ (#|^, L (0)|vac). (19) 

This condition leaves one phase ambiguity: there is still an invariance under the change 
D±(9) i — ^ u ±1 D±(9) for \u\ = 1. This phase ambiguity will not play any role in the following, 
hence we fix it arbitrarily by choosing the solution (12), (14). 

2.2 Primary twist fields 

The Dirac theory enjoys a U(l) internal symmetry: ^ r.l e 27Tla ^R.L, oi S [0,1). The 
associated primary (with respect to the Fermion algebra) twist fields, a a (x,y), are local, 
spin-less, U(l) neutral quantum fields with bosonic statistics, and are known to have scaling 
dimensions a 2 [20]. Their twist property, characterised by a locality index e 2ma , can be 
expressed via a monodromy property of correlation functions of time-ordered operators, 

C{z) := (vac|T[---*fl,i(ar,3/)£7 a (0)---] |vac). 

Seeing C(z) as a function of the complex variable z smoothly continued through the y = 
line, its continuation along a loop 7 surrounding once counter-clockwise the origin is given 

by 

C{e 2m z) = e- 27ria C(z). (20) 

This is true for any loop 7 that can be contracted to the origin without intersecting the 
position of other twist fields possibly present in the correlation function. A similar relation 
holds for the Hermitian conjugates "f^, L , but with the factor e 2ma on the right-hand side, 
instead of e~ 2ma . This twist property can also be expressed operatorially via equal-time 
exchange relations: 

f (7 a (0)^fR, L (x) (X<0) 

* r L(x)cr a (0) = < (21a) 
f «7a(0)¥i,(aO (a? < 0) 

*U(*K(o) = <^ ■ \ (21b) 

1 e- 2 -M0)*U*) (*>0). 
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These naturally lead us to define primary twist fields with negative index via Hcrmitian 
conjugates: 

e[0,l). (22) 



These exchange relations along with the property of being primary immediately gives the 
form factors of the U(l) twist fields [11, 16, 8]. As the fields are U(l) neutral they only have 
non-vanishing form factors with U(l) neutral states, so for —1 < a < 1, 

(vm\a a (0)\e 1 e 2 ---9 2n ) + ... + -...- = c a m°? {^^f (ll2=i(«i)* +a («i+«)*- a ) 



n"=i Wii n+ i(u r +u s 

where there are n — 's and n +'s in the state, and Ui := cxp(#j). In particular, the two-particle 
form factor is 

2sinf7ra) p a ( 9l_6j ) 
(vacM0)|M 2 >+- = c a m a ^ C och e ^ " ( 24 ) 

Note that the form factors (23) may differ from those written in other publications in the 
choice of the sign of a and in the choice of the sign of the two-particle form factor (the 
latter amounting to an overall sign (—1)™). These choices relate to our choice of exchange 
relations (21) (defining the operator er Q ), and to our choice of crossing symmetry relation 
(19) (defining the asymptotic states). Calculations below, using the two-particle form factor, 
will show that the two-particle form factor (24) leads to a behaviour of ^ R (x)a a (0) and 
\& n(x)(T a (0) as x — ^ exhibiting a branch point in accordance to the exchange relations 
(21). Further calculations will show that the disconnected, delta- function term in a crossing 
symmetry relation similar to (19) is correct, this being rue only with the sign chosen for the 
two-particle form factor. 

The constant c a is chosen, up to a phase, in order to guarantee the CFT normalisation 
of the fields: 

(vac|cr a (a;i,yi)(T_ a (a;2,y2)|vac) ~ 2 as \z\ — z?\ -t Q. (25) 

\zi - z 2 \ 2a 

The phase of c a is fixed by requiring that it be real and positive. The evaluation of the 
constant c a was first explained in [15] (this constant was first obtained in an unpublished 
note by ALB. Zamolodchikov). This constant, or rather its generalisation to the Dirac Fermi 
field on the Poinare disk, was identified with Barnes' G- function in [3]: 

Ca = G(l-a)G(l + aY (26) 

Recall that the main property of Barnes' G-function is G(l + z) = T(z)G(z), with the 
normalisation G(l) = 1. In the next section we provide an alternative way of evaluating c a . 

Other matrix elements can be obtained by analytic continuation in the rapidities, as is 
standard in the context of 1+1-dimensional QFT [23]. Note that the form factors (23) and 
the Hermiticity relation (22) are in agreement with the analytic-continuation formula 

(vacMO)|(0i + wr) • ■ ■ (0 2n + wr)> + ... + _..._ = +-+"—<0 an • ■ ■ 0iK(O)| vac), (27) 

where the analytic continuation is simultaneous on all rapidities. 



(23) 
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3 Descendent twist fields and their form factors 



We are interested in deriving differential equations for two-point correlation functions of 
the primary twist fields a a . However, in the derivation, the intermediate steps naturally 
involve other twist fields; descendents with respect to the Fermion algebra. Here we provide 
unambiguous definitions for these fields and evaluate their form factors. 



3.1 Basic definitions 

First, notice that (23) is an analytic function of a on C \ Z*, for fixed rapidities; the points 
Z* := Z \ {0} correspond in general to poles. Hence, it is natural to define the fields <r a 
for all a £ M. \ Z* by formula (23) (and by analytic continuation in the rapidities for other 
matrix elements). Then, the two-point function of such twist fields is, for any fixed distance 
between the fields, an analytic function of a on a neighbourhood of M. \ Z*. In fact, we 
expect that the coefficients of operator product expansions with the stress-energy tensor are 
"uniformly" analytic in a at all distances. This means that the new field defined by analytic 
continuation has dimension a 2 ; we also expect the correct CFT normalisation (25). 

Other descendent twist fields, of Fcrmionic statistics (in particular, whose form factors 
with even number of particles are zero), can be obtained as regularised limits where the 
position of a Dirac field approaches that of a twist field a a ; equivalently as coefficients in the 
operator product expansions (OPEs). For a in particular subsets of R \ Z*, we define two 
one-parameter families of Fermionic twist fields via the coefficients occurring in the leading 
short-distance asymptotics associated to the Dirac fields and "J^,: 

T[& R (x, y)a a (0)] ~ (-iz) a a a+ha (0) for a < l - 

T[^ R (x,y)a a (0)}r,(-iz)- a a a . 1 , a (0) for a > ~. (28) 

As above, T is the time-ordering operation. The power functions on the right-hand sides are 
on their principal branches; note that there is agreement between the phases occurring from 
exchange relations (21) and the phase differences occurring through the cuts of the principal 
branches. The requirement of phase agreement is not quite enough to fully determine the 
particular power functions occurring: these are consequences of CFT considerations (see 
appendix A), or of form- factor calculations (performed below). From (28), we see that the 
new fields a a ±\^ a have dimensions a 2 ± a + 1/2, spins ±a + ^ an d charges =pl. They also 
satisfy the hermiticity relation 

a a ±i,a = c-c^l >-«• ( 29 ) 

The range of a shown is a consequence of the equations of motion, as we explain shortly. We 
define fields for all a £ R \ Z* again by analytic continuation beyond the range of a shown. 

CFT considerations (appendix A) also suggest that we may define the same Fermionic 
twist field families using and ^>' L : 



T 



y{(x,y)a a (0) ~-(iz) a cr Q , ia _i(0) for a > ~ 



T [*£(*, tfK(0)] - -(iz) a a a , a+1 (0) for a < ^. (30) 

Together with the equations of motion (9), this suggests the form of the leading OPE for all 
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a£R\Z*: 



T 



(-iz) a a a+1 , a {0) + 



and 



T[*«(x,y)£r a (0)] 



T [^(x,2/)a Q (0) 
T[*i,(a;,2/)(r Q (0)] 



(-iz)- Q a Q _ ljQ (0) + 

-(w)- Q ^,a-l(0) 

-(iz) a a a . a+1 (0) - 



1 - a 
in 



1 + a 



{izy- a a a , a -i{0) 
{iz) 1+a a aia+1 {0) 



(31) 



■(-«) 



' a+l,a 



(0) 



1 + a 

{-iz) 1 -^-!^), (32) 



1 — a 



where the meaning of ~ is that for any value of a, the term on the right-hand side that is 
dominant gives the leading OPE. Note that according to these equations, the leading term 
of the OPE occurring in the CFT is modified in certain ranges of a by lower-dimensional 
"ghost fields", like ma ai0l —i, which make no contribution in the massless limit. Note also 
that sub-leading terms in these OPE's will have poles at other integer values of a. 

Moreover, further applications of Dirac fields may lead back to the original family a a of 
twist fields. For instance, we have, for a < 0, 



T 



$ R (x, y)a a , a +i (0) ~ -i {-iz) a cr a+1 (0). 



(33) 



This relation is an easy consequence of CFT considerations (see appendix A). Note in par- 
ticular that the double application of Dirac fields preserved the CFT normalisation, because 
the Dirac fields themselves are CFT normalised. 

The next two subsections are devoted to computing the form factors of the Fermionic 
twist fields defined in (28), and to showing the correctness of (31) and (32). In combination 
with these form factors, the consistency of (28) and (30) leads to the non-trivial recursion 
relation (49) for the normalisation constant c Q . Equivalcntly, relation (33) can be used 
(since (33) is a simpler consequence of CFT, this is a slightly more direct way of obtaining 
the recursion relation (49)). We will show that it is indeed satisfied by (26). Interestingly, 
such a recursion relation can be seen as a novel way of evaluating this constant. 



3.2 One particle form factors 

In this subsection we give details of the computation of (vac|<r Q +i iQ (O)|0) + . The form factor 
is evaluated by expanding (va,c\^> R (x)<r a (O)|0)+ and finding the term proportional to z a , 
using (28). We provide the result for the form factor (vac|cr Q _i jQ (0)|6')_, which was obtained 
by expanding (va,c\ 1 ^ft,(x)a a (0)\8)^ similarly. 

In order to obtain the leading short-distance behaviour of 

(vac\tf R (x,y)a a (0)\6) + (34) 

we evaluate it by inserting (18) between the two fields. We consider x < 0, y — for 
simplicity, as this is sufficient to extract the onc-particle form factors. The full z dependence 
of the leading behaviour can be restored using rotation (or relativistic boost) covariancc. 
There is only one non-zero form factor of ^ R given by 

(v&c\y R {x)\0) + = V^e e/2 e lxpe (35) 
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This means that there is only one term involving a a to be evaluated, which can be accom- 
plished using the form factor rules of [23] and the two- particle form factor (24): 

+ (</>|cr Q (O)|0) + = (vac|er a (O)|0,</> + w-iO+) + _ + (vac|cr a (0)|vac)<S((9 - 0) 

= Cam a2 !^ e — e ?°~l 0+ s + c a m a2 S(e - 0). (36) 

2n sinh ( e ^+' 0+ 

The +i0 + prescription is to be understood in terms of distributions: it is a prescription as 
to how an integral should avoid the pole of the function. Putting these together we see that 



2 

+ c Q m a2+1 /V/V^. (37) 

The integral is conditionally convergent. In order to make it a convergent integral, we make 
a change of contour, shifting (f> t— > <fi — in/2. This shift of contour does not give any pole 
contribution thanks to the i0 + prescription, and the result is indeed a convergent integral 
for x < 0: 

(vac|*kMK(O)|0> + = Cam ^i/^ e - m /2| d ^-./4 e tf . ff' 

siring 2 + 4 J 

+ c a m a2+1 / 2 e e / 2 e ixpi> 

i e 2 — e 2 

+ c a m a2+1 / 2 e 9 ' 2 e lx P<>. (38) 



In general, the leading behaviour at small x will be obtained from the large-</) behaviour 
of the integrand. Let us analyse what happen at x = in different regions of a. In the 
regions a < and a > 1, the resulting integral is divergent. In these cases, the leading 
small- a; behaviour is growing, and can be obtained from the leading large- \<f>\ behaviour of 
the integrand only. On the other hand, in the region < a < 1, the integral converges, and 
the resulting constant is cancelled by the term arising from the delta function (see Subsection 
3.3). Note that this cancellation requires in particular the correct sign of the two-particle 
form factor (24). Hence, in this case, the leading small-x behaviour is decaying, and can be 
extracted again from the leading large- \4>\ behaviour by first rendering the integral divergent 
by taking its first ^-derivative. As explained before, we expect the field cr Q+1 a to occur in 
the leading OPE for a < 1/2, and to have an appropriate ghost fields in the range a > 1/2, 
as per (31). We will see below that the one-particle form factors agree with this. In order 
to evaluate the form factor, it is sufficient to consider the region a < 0. 

In order to extract the small- a; behaviour in the region a < 0, we split the integral into 
two regions, </> € (— oo, 0) and <fi £ (0, oo). In each region we pull a factor out of the fraction 
to leave a denominator of the form 1 + r/ where |ry| — > as \<p\ — > oo, and we then Taylor 
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expand each fraction. This gives 
(vac|*Ux,O)<7 Q (O)|0) + ~ c Q m Q2+1 / 2 ^^ e -^/ 2 x 

dcj, (_ e (a+l/2)0- a <p+xE 4l _ ie (a+3/2)8-( a +l)<P+xE 4 , + ^ _|_ 



_ ie (a-l/2)e-{a-l)<P+xE <l> _ e (a-3/2)8-{a-2)j>+xE< l , 



(39) 



The leading divergent term for a < occurs in the region (0, oo). This term can be rewritten 
using a modified Bcssel function: 

-atj}+xE,p\ I Jl f ~ a4>+mx cosh 4>\ _ I j ± I — aip+mx cosh <j>\ 



d(j) (e~°"*' +x *) = / d</> ( e - Q< P+ m:rcosh< P) — / d0 (g-acp+mxcoshi^ 
J —oo J — oo 

= 2K\ a \(-mx) + 0(1) (40) 
as x — >• 0~. Hence, the latter integral diverges like 

r(-a)( — ) °. (41) 
\ —mx J 

The first equation in (28) allows us to identify the coefficient of the divergence with the 
one-particle form factor we are seeking. We obtain: 

(vac|cr a+1)a (O)|0> + = Ca i-!!^ m " 2 +"+V2 e (a+l/2)^ (42) 
T(l + a) 

By analytic continuation, this provides the form factors for all a € R \ Z* . 

In order to verify the full z-dcpcndcnce of the leading behaviour (28), we may use the 
spins of the operators involved and examine the effects of a relativistic boost of rapidity (3 
about the origin. Acting with the boost operator Rp, defined essentially by Rg\9) = \0 + /3), 
we see that 

Rptf R {x,y)Rp = e^H R {xp,yp) and Rpa Q (0)R} 3 = a a (0) (43) 

as the fields have spins 1/2 and respectively. Here, xp and yp are obtained using zp = e^z 
and zp = e~"z. Moreover, from the form factor (42), we find that 

Rpa a+1 , a (0)Rl = e<- a+1 /^a a (0). (44) 

Hence, the field cr Q+ i, Q has spin a+1/2 (as mentioned above). When acting with the rotation 
operators on ^ R {x, j/)<r aja (0), we can either perform the boost first and then use the OPE: 

Rp^ R (x,y)a a (0)Rl = e fi/2 ^ R (xp,yp)a a (0) 

r,e^ 2 (-ize^) a a a+1 , a (0) (45) 



or vice versa: 



Rp^ R {x,y)a a (0)Rl ~ {-iz) a Rpa a+Xja {Q)R\ 

= (-iz) a e<- a+1 ^a a+ha (0)- (46) 
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Agreement confirms that the first equation of (28) is satisfied. 

The one particle form factor of the field er Q -i, Q (0) can be calculated in a similar way, 
using the second relation of (28). The calculation, in the range a > 0, gives 

(vac|a Q _ liQ (0)|#)_ = -ic*4\ - m « a -«+ 1 /^- a+1 ^ 6 (47) 

r(l-a) 

which, by analytic continuation, holds for all a 6 ffi \ Z* . Note that 

(vac|cr a _i ia (O)|0)_ = -«(vac|cr_ a+ i 1 _ a (O)|0} + . 

With the crossing symmetry relation (vac|<7a._i iQ .(O)|0 + Z7r)_ = + (#|cr Q _i, Q (0)|vac) along 
with the hermiticity relation (29), this implies 

-i(vac|o- a ,_i )O ,(O)|0 + wr) + = ((vac|cj a _i ia (O)|0) + )* 

which is indeed satisfied by (42). 



3.3 A recursion relation for c a , and further verifications of short- 
distance behaviours 

We now proceed to verify the relations (31) and (32), as well as (33) and other similar 
relations. As we will show, these in fact provide a non-trivial recursion relation for the 
constant c Q . This recursion relation can in fact be used to evaluate the constant c a , since, 
for instance, relations (31) were deduced solely from CFT arguments. We note, in particular, 
that this way of evaluating c a implements the CFT normalisation of the twist fields somewhat 
more clearly than in previous methods [15]. 

First, we may verify the first relation of (31) in the range a > 1 by evaluating the leading 
diverging term of (va,c\^\ i (x)a a (0)\8) + for a in that range. Starting with (39), we need only 
to consider the region of integration (-co, 0). We can re- write the integral of the first term 
displayed using a modified Bessel function, in a similar fashion to what was done in the 
previous subsection. Comparing with the mass term in the first relation of (31), we find the 
equation 

— iira/2 

(vac|a Q . Q _ 1 (O)|0) + = ic a ^- -m a — +V2 e «(a-i/a) (48) 

1 (1 — a) 

After a shift a h-> a + 1, comparison with the form factor (42) gives the relation 

ja_ = r^ + i) 

c a +i 1 (-a) 

It is a simple matter to verify that this relation is indeed satisfied by (26) using the properties 
of the Gamma function and of Barnes' G-function. We remark that this provides a very non- 
trivial check of the validity of the definition of the fields u a ±i,a by analytic continuation. 

The recursion relation (49) along with the symmetry property c a = c_ a gives a way of 
evaluating c a . Indeed, although these relations have many solutions, the solution is made 
unique by imposing additionally the normalisation condition cq = 1 and the conditions of 
monotonicity and convexity in the range a € [0, 1). These extra conditions are indeed 
satisfied by (26). The normalisation condition Co = 1 is a simple consequence of the fact 
that (Tq is the identity operator at a = 0. However, we do not yet have clear arguments for 
monotonicity and convexity from quantum field theory. We note that all calculations in later 
sections are unaffected by the explicit choice of the solution c a , and only follow from the 
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recursion relation (49) (along with, implicitly, the symmetry property and the normalisation 
condition). 

Second, we may complete the verification of the first relation of (3f) by analysing the 
region < a < 1 . We first note that the x = value of the integral in the first line of (38) is 

„a(S — <p) n ma/2 

e = _^^ e */2. (50) 

sinh(0_^ + Hl) sm(Tra) 

This can be obtained by first shifting <fi i— > <f> + 9 in order to extract the 0-dependent factor, 
then by shifting the contour cf> i-> <fi — 2in, getting the pole at <f> — — 3i7r/2, and noticing 
that the integral with a shifted contour is just the initial integral up to a phase (solving the 
resulting equation gives the answer). Hence, we indeed find that the quantity (38) decays 
as x — > 0. The power law of the decay can be obtained simply by taking a derivative with 
respect to x: it is then a divergent behaviour, hence can be obtained by the techniques 
illustrated above, using (39). We then integrate over x in order to recover the leading OPE. 
The leading term is different in the regions a < 1/2 and a > 1/2. For a < 1/2, we recover 
the first equation of (28), with (42). For a > 1/2, taking into account the conditions on 
the spin, the power law is (iz) 1_a , and the coefficient is identified with a multiple of the 
one-particle form factor of the field ma a , a -i- The proportionality constant can be evaluated 
using (26), and this again only involves verifying the recursion relation (49). The resulting 
leading OPE is in agreement with the first equation of (31) in the range 1/2 < a < 1. The 
case a = 1/2 simply involves the two terms calculated, again in agreement with (31). 

Third, in a similar way, we may verify agreement of the form factor (47) with the second 
relation of (31) by considering the ranges a < — 1 and — 1 < a < 0, and by using similar 
arguments as those above. For the range a < — 1, the mass term of the second relation of 
(31) leads to the equation 

p iwa/2 

(vac|a a , Q+1 (O)|0)_ = c a —— —m a +«+i/2 e *(-a-l/2)_ (51) 
T(l + a) 

Again, comparing with (47), this leads to the recursion relation (49), hence is satisfied by 
(26). Likewise, the relations (32) can be verified along entirely similar lines. 

An alternative way of deriving a recursion relation for the constant c a is to notice that 
multiple applications of the Fermi fields can lead back to the family of primary twist fields 
<T a . By rearranging (33) we have that 

(T Q+ i(0) - i{-izy a T[^ R {x, j/K, Q+ i(0)] (52) 

as z — > (for a < 0). The vacuum expectation value of cr a +i(0) can be calculated from 
this relation in the same way as the one particle form factors were calculated in Section 3.2. 
Firstly (18) is inserted between the fields to give 

(vac|^(x,y)cr a , a+ i(0)|vac) = J d9 (vac|*| ? (x,?/)|6l) + + (6l|CT Q , Q+ i(0)|vac) 

=Vm J d9 e <>/2 e -vEe+ixpe (vac|CTa,a+l(O)|0 + in) - 

= - ic a m a2+a+1 C , " a/ \ [ d9 e- ae - yEe+txpe . (53) 
T(l + a) J 

Doing as before, we set y = and shift the integration contour, 9^9 — iir/2. We then take 
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the small x expansion to obtain: 



(vac|*kMK, a+ i(0)|vac) = -ic a m a * +a+1 — j d9e- a9+xE » 

_ ;„ _a 2 +g+i T (- a ) ( 2 



r(i + Q) 

where on the second line we assumed a < 0. Hence (52) gives 



(54) 



(vac| ( r Q+1 (0)|vac) = c Q+1 m (Q+1)2 = c a m a2+2a+1 J} a) (55) 

T(a + 1) 

which gives again the recursion relation (49). 

Finally, it is worth mentioning that the matrix elements (v&c\a a (0)^ R L (x)\6) e and 
(vac|fj a (0)\I/ r,l{x)\0) € can be evaluated by similar methods, but without the need for us- 
ing crossing symmetry. The result is in agreement with the exchange relations: for in- 
stance, for x < 0, the same function of x and 9 is obtained as for the matrix elements 
(vac|^I / ^, l (x)<Tq,(O)|0) £ and {v&c\$! r^l{x)&u (O)|0) c respectively. 



3.4 Higher-particle form factors 

All higher-particle form factors can be obtained by Wick theorem. The form factor 

(vac|o- a± i iQ (0)|6>i, . . . ,6>2„+l)ei,...,e 2n+ i 

is the sum over every term obtained by a product of n (normalised) two-particle form factors 
of a a associated to n disjoint pairs of particles (contractions), multiplied by the one-particle 
form factor of cr a ±i ta associated to the remaining unpaired particle, with an appropriate sign 
depending on the number of crossings of the contractions: 

(vac|a- Q ±i jQ (O)|0i, . . ■ ,6 , 2n+i)e 1 ,...,£ 2 „ +1 (56) 
£ ( _ ir o„ Yl ' al ;'jJ 2 tj ^ (y^\a a± U0)\0 q )e q . 

partition {q; (sj,tj):j=l, ...,»} VJ^ 1 " / 

of {!,... ,2n-l} 



4 The double model 

In order to make progress in calculating two-point correlation functions of the twist fields wc 
take inspiration from [9] where two non-interacting copies of Ising field theory arc considered. 
Since each copy is quadratic, there is invariance under rotations amongst the two copies, 
which is an extra continuous symmetry. In the Ising model, the associated conserved charge 
allowed differential equations for the two-point spin-spin correlation function (and for other 
objects) to be written down. 

We consider two copies of the Dirac model defined in Section 2, which we denote by 
^ and $. These copies have the creation and annihilation operators D±\6) and E±\o) 
respectively, and independently satisfy all the relations of Section 2. The two copies do not 
interact, hence obey the anti-commutation relations: 

{*% l (z),* 1 £lW}=0 and {D ( 2\9i),E ( i ) (6 2 )}=0. (57) 
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There are also independent twist fields associated to the two independent U(l) symmetries 
of the two copies. These twist fields will, when relevant, carry an extra superscript to denote 
the copy they belong to: for instance, <r* and cr*. When convenient, we will use such 
superscripts more generally to denote operators belonging to the respective copies. 

4.1 Conserved charges 

Energy and momentum are conserved charges of the single-copy model. They are associated 
to the dynamical invariance under space and time translation: the equations of motion (9) 
possesses this invariance. In the double-copy model, the energy and momentum conserved 
charges are the sums of the corresponding charges of each copy. However, there are other 
conserved charges in the double-copy model that can be constructed out of these. As the two 
copies are non-interacting, the energy-momentum operators for each are still independently 
conserved quantities, and can be combined differently to give new conserved quantities. The 
two specific conserved charges that are of interest to us are the differences of those of the 
single copies. We define P and P via the following action (these conserved charges are chosen 
to be anti-Hcrmitian) : 

[P, 0*0*] =idO*Q* - lO^dO* 

[P } o*o*] = _iao*o* + io*ao*. (58) 

This holds for any local fields 0* and 0* interacting non-trivially with fields in copy and 
$ respectively In particular, the actions on the creation and annihilation operators are 

[P,D±{0)} = -ime e D±(6), [P,E±(6)] = ime e E±{9), 
[P,D±{6)\ = ime- e D±{6), [P,E±{6)\ = -ime~ 6 E^{0) 

from which it is simple to derive an explicit expression for P and P through bilinears in the 
creation and annihilation operators. 

In the double-copy model, there is yet another conserved charge Z, related to the 0(2) 
rotation symmetry amongst the copies. Written in terms of the Fermi fields W and $, the 
charge Z is (again, chosen to be anti-Hcrmitian) 

Z = h\ Ax{ ^ R ^R + * L$ i + + (59) 

The action of this charge on the creation and annihilation operators is 

[Z,Dl(6)] = -Et(0) [Z,El(9)]=Di(9) (60) 

and similarly for D±(8) and E±(8). 

It is of course possible to obtain higher order conserved charges by commutations amongst 
P, P and Z. In the approach of [9], where the Ising model was studied, it is assumed that 
P, P and Z all preserve the quantum state where correlation functions are evaluated, and 
the Ward identities corresponding to certain of these higher order conserved charges are 
used to derive differential equations for twist-field correlation functions (it was also observed 

there that the higher order conserved charges form a si (2) algebra). However, it turns out 
that this is not necessary, as was noticed in [4]. Indeed, the Ward identities coming from 
the invariance under the Z-action alone are sufficient to obtain non-trivial equations for 
correlation functions - we do not need to assume that the dynamical space-time translation 
symmetry is a symmetry of the quantum state under study. When the quantum state does 
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possess this symmetry, the space-time dependence of correlation functions simplifies, and the 
resulting equations are the integrable differential equations describing the Ising correlation 
functions. As we will see, the same structure occurs in the present case of the Dirac theory 
(with the additional condition that we need parity invariance in order to further simplify the 
space-time dependence of correlation functions). This is an important remark, as it shows 
that the method may have applicability beyond the cases where the quantum states are 
invariant under space-time translation (and parity). 

Below, we will derive the integrable differential equations from the Ward identities asso- 
ciated to Z, and use the simpler space-time dependence of correlation functions coming from 
translation and parity symmetries. We refer to Appendix B for the equations resulting from 
the Z Ward identities when these symmetries are not present. 

As was observed in [4] (in the generalised situation of a theory on the Poinarc disk) , it is 
the relation 

[P,[P,Z]] = [P,[P,Z]]=Am 2 Z (61) 

which provides the dependence on the mass in the differential equations. Relation (61) 
corresponds essentially to the equations of motion of the theory, as expressed using the 
charge Z. It can be derived from (58) and (59) using the equations of motion (9). Relation 
(61) is the only one where the massive theory is used: besides it, we only need the action 
of Z on products of twist fields and their first derivatives in order to find the differential 
equations, and this action, as is verified below, is solely a consequence of CFT calculations 
(i.e. calculations using the massless limit of the theory). 



4.2 Equations for the action of Z 

The action of the conserved charge Z on products of local fields 0*0* gives linear com- 
binations of similar products of local fields (with unchanged locality index) whenever the 
combined locality index of the product 0*0* is one. Otherwise, the result of the action 
is non-local, and the resulting Ward identities arc not useful. We expect that the relations 
derived in the CFT limit in Appendix A, where the right- and left-moving components of 
the Fermions are separated, will carry across to the massive case when written in terms of 
well-defined twist fields of the massive theory. Although we have not provided a full analysis 
of this fact, which would involve generalising the OPEs derived in Appendix A to massive- 
theory OPEs involving mass terms as in (31) and (32), we have nevertheless verified it using 
form factors with up to 4 particles. 

The equations representing the action of the charge Z on twist fields involving no deriva- 
tives are: 

[2,<M = (62a) 

[Z,(7 a a^_ 1 ] = i(<7a_i )ai er a , a _i - ^Z,a.-l a a-l,a) ( 62b ) 
[Z, cr*cr* +1 ] = i(cr* + i ia cr* Q+1 - (T* Q+1 cr* +l Q ) (62c) 
[Z,o-* +liQ er* Q+1 ] = i(er*er* +1 - er* +1 cr*) (62d) 

[ Z i a t+l a t-l\ = — (( 9f7 *a+l) CT a,a-l ~ a t,a+l( da t, a -l) + 

(^a+l,a) CT o-l,a - CT *+l,a(<9CT*-l,J) (62e) 
[Z^t-l^t+l] = - (^a,o-l( 9o '*o+l) - ( 9<7 a,o-l) CT o,a+l+ 

°"*-l,a(^*+l,a) - (^a-l,a)°"J+l,a) ■ ( 62f ) 
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Here and in the rest of this section, all fields arc evaluated at one single point, say for instance 
the point (0,0). When verifying these equations we place the vacuum state on the left and 
a multi-particle state on the right. The left-hand side is expanded using (60) along with 
(vac|Z = and Z\va,c) = 0, this gives form factors of the type (23). The right-hand side 
gives form factors of the type found in Sections 3.2 and 3.4. The two sides are found to be 
equal. We note that for this equality to occur, we do not need explicitly the form (26) of 
the vacuum expectation values, but only the recursion relation (49). This agreement is very 
non-trivial and is strong evidence that Equations (62) hold for all matrix elements. 

The origins of Equations (62) are not obvious without reference to the CFT. Equation 
(62a), however, is easy to verify, and using this (62b) and (62c) can be 'derived' from the 
relations of the type 

*fl*LCT Q i-> (T Q+ i and ^ B ^ L a a i-> <r Q _i 

(see (30) and (33) for the more accurate relations). The final three equations do not appear 
to have simple origins outside of the CFT but can be verified using the method described 
above. It is also worth noting that the left-hand sides of (62b) and (62c) are the same up 
to a shift in indices and a relabelling of the fields. The right-hand sides are also consistent 
with these operations. 

Equations involving one derivative are as follows, and have been derived from CFT and 
verified using form factors like those above: 

[Z, (da a )cr a \ = a(<J a _ l a <7 a+l a -<T a+l a a a _ ha ) 
[Z, [da a )a a \ = a{a a a _ 1 a a a+1 -a a a+1 a a a _ 1 ) 

[ Z , (^aVa+l] = '((^a+O^+l.a ~ ^l.J^a+l)) 
[ZAdVa^t+l] = ^a+l^a+l.J ~ (^*+l,a) CT *a+l) 
[Z, {dat)°t-l\ = i (( 5CT *a-l) £r J-l,a- £r a-l,a( 9o '*a-l)) 

[s,(&7>2-i] = >X< a -i(^i,J-(^i,J< a _ 1 ) 

[ Z , K,a+>!+l, J = K a a+l( 9ff t) - ( 9cr a Va+l) 

[Z, (9o-*+l, a )o-* a +l] = K a a( d ^+l) ~ (^a+lVS) 

[Z, (^U+lVj+l.d = KidVa+l)^ -0"*(9<7* +1 )) 

(9o-* +1)Q )<7* a+1 ] = i((9cr*)cr* +1 -CT* +1 (<9cr*)). 

Note that in Section 5, we only need few of the zero- and one-derivative equations showed 
here. 

Finally, the action of Z on fields involving higher numbers of derivatives are not in 
general expected, in the massive theory, to be in agreement with calculations from CFT. 
However, the only ones that we actually need can be obtained by using those above along 
with the equation of motion (61) of the charge Z. Indeed, we only need actions of Z on 
double-derivative fields of the form [P, [JP, 0*0*]], which we can evaluate using 

[Z, [P, [P, 0*0*]]] 

= [[P, [P, Z}], 0*0*] + [P, [Z, [P, 0*0*]]] + [P, [Z, [P, 0*0*]]] - [P, [P, [Z, 0*0*]]] 
= 4m 2 [Z, 0*0*] + [P, [Z, [P, 0*0*]]] + [P, [Z, [P, 0*0*]]] - [P, [P, [Z, 0*0*]]] 

This is how the mass dependence will appear in the equations for correlation functions. 
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5 Correlation functions 



5.1 Correlation functions of interest 

The aim of this section is to derive differential equations for two-point correlation functions 
of twist fields using Ward identities associated to the conservation of the charge Z. We will 
keep the state where the correlation functions is evaluated as general as possible. 

The first general form of a state invariant under Z is the mixed state where the average 
of an operator A is given by 

tW' (63) 

where V is any bilinear form in the creation and annihilation operators (or in the Dirac 
Fermi fields). In the double copy model, it is clear that + is invariant under 0(2) 
rotations, hence that 

Tr (e v *+ v *A) 



(64) 



is zero whenever A is of the form [Z, ■ ■ ■ ]: these are the Ward identities associated to Z. These 
states include, for instance, the finite-temperature state V = —H/T, with temperature T, 
as well as the vacuum state, the limit T — > 0. 

Another form of the quantum state is that where averages are given by 

(B'\A\B) (65) 

where \B) and \B') are boundary states of the form \B) = e y |vac) and \B') = e v |vac), for 
bilinear forms V and V in the creation and annihilation operators. Again, 

((B'fiB'flAlB^B*) (66) 

is zero whenever A is of the form [Z, ■ ■ ■ ]. Such states include, for instance, the integrable 
boundary states [10] of the Dirac theory. 

Further, since the action of the Z operator on local fields is a local property, it is not 
affected by the exact space where the quantisation holds: we may also choose a Hilbert space 
of quantisation on a finite segment (e.g. the circle of length L), and again consider the states 
mentioned above. 

We will use the single notation (A) in order to represent averages (63), (64), (65) or (66). 
We will show that the sinh-Gordon partial differential equation determines the following 
three families of correlation functions: 

F(x,y) = (a a (x,y)a a (0)), 
G(x,y) = (a a (x,y)a a+1 (0)), 

H(x,y) = (cr a+1 . a (x,y)a a . a+ i(0)) (67) 

(68) 

whenever the correlation functions are translation and parity invariant. Here, the dependence 
on a is implicit, in order not to clutter the equations. In the same spirit, we will also use 
the notation 

F(x,y) = (a a+1 (x,y)a a+1 (0))). (69) 

When there is no translation or parity symmetry, more general equations can be derived, see 
Appendix B. 
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We note that correlation functions involving different ordering of the two fields present 
in the definitions of F(x,y), G(x,y) and H(x,y) are simply related to those written above. 
This is obvious for the cases of F(x, y) and G(x,y), since the twist fields involved commute. 
For the case of H(x, y), we must apply the braiding relations (21). Recall that by definition 

a a+ i a (x,y) = lim (-i(z - w))~ a ^ R (x' ,y')a a (x,y) 

a a , a+ i{x,y) = ]im(i(z- w)y a ^ L (x' ,y')a a (x,y) (70) 

z — >w 

where z = — ^(x' + iy') and w — ~^(x + iy) are appropriately time ordered (y' > y). 
The factors (—i(z — w))~ a and (i(z — u>))~ Q arc taken on the principal branch, and so are 
continuous exactly where ^ R (x', y')a a (x, y) and ^l(x', y')a a (x, y) are continuous. Using 
these definitions and (21) we obtain 

^p,p+i{X2,y2)(T a+ i, a {x 1 ,y 1 ) = < _ e2m if < >a a+ i, a {x 1 ,y 1 )ap^ +1 (x2,y2)- (71) 



5.2 Translation and parity symmetry 

Space and time translation invariance simplifies the functional form of the correlation func- 
tions. It implies that it is indeed sufficient to fix the position of the second field to in all 
correlation functions above, as we have done, and further that 

F{-x,-y)=F(x,y), F(-x,-y) = F{x,y). (72) 

It also implies, using the fact that the fields a a and cr Q +i commute and using the braiding 
(71), that 

(<T a+1 (x,y)a a (0)) = G(-x,-y), (a ata+ i(x,y)a a+ltCt (0)) = -e 2m H{-x,-y). (73) 

Space and time parity symmetry (i.e. invariance under simultaneous inversion of space 
and Euclidean time) further simplifies the functions G and H. In two dimensions, this is 
equivalent to Euclidean rotation by an angle of tt. Since this is a Euclidean rotation, some 
subtleties occur. The symmetry cannot be expressed simply through a unitary operator, 
but rather is a relation connecting vacuum expectation values. It says that there exists a 
unitary operator R, which performs a parity symmetry on individual local fields, such that 
the following equality holds: 

(vac|Oi(o;i,yi)---0„(a;n,yn)|vac) = (vac|i?C„(x„, ?/„)■■■ Oi(xi, yi)R^\va,c) (74) 

whenever the local fields Oj(xj 7 yj) are time-ordered: y± > ... > y n . The unitary operator 
does not act trivially on the vacuum, because it transforms an in-vacuum into an out- 
vacuum and vice versa. But it acts locally on local fields, in particular inverting both space 
and time coordinates, so that the correlation function on the right-hand side is still time- 
ordered. Relation (74) has its origin in the Euclidean path-integral picture, whose connection 
to the operatorial picture requires time-ordering: it represents the invariance of Euclidean 
correlation functions under a rotation of the Euclidean plane by tt. 

The operator R can be written explicitly. Naturally, it could be seen as a shift of rapidity 
by the imaginary number in. However, there is another, more appropriate way of representing 
this operator, with the same effect on Fermion fields: the operator R can be seen as a 
combination of two symmetries, space parity S and time parity T, with 

R = TS. 
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Space parity is implemented by the unitary linear operator S satisfying 

S* R>L (x, y)S = ±M L , R (-x, y), S 2 = 1, (75) 

and time parity is implemented by the unitary linear operator T satisfying 

T9 RiL (x,y)T = * L , R (x,-y), T 2 = 1 (76) 

(note that the latter is not an anti-linear operator, because it inverts the Euclidean time 
instead of the real time - it is akin to a parity operator). It is easy to verify that these are 
indeed dynamical symmetries: they preserve the equations of motion (9). The actions of S 
and T on the mode operators D±(6) can be obtained from the relations above, and are given 
by SD±(6)S = TD±{-9) and TD±(9)T = -iD^(-6). The combined parity transformation 
is 

TS^ RiL (x,y)ST = ±i^ RiL (-x,-y), TSD ± {9)ST = ±iD%(9), (77) 
and it has the following simple operator representation: 

(78) 

The transformation properties (77) tell us that the creation operators are transformed 
into annihilation operators of opposite charges, and vice versa. In combination with (74), 
this allows us to extend the invariancc property to correlation functions in arbitrary states: 
we have in general, again with time ordering y\ > . . . > y n , 

{B'ldfa, yi) ■ ■ ■ O n (x n ,y n )\B) = (B\RO n (x n , y n ) ■ ■ ■ O^x^y^R^B') (79) 

where \B) and \B') are the states obtained from \B) and \B'), respectively, by inverting the 
charges of all particles, with extra factors —i and +i for each particle if the new charge is — 
or + respectively. The action of a rotation by 7r on multi-particle states is indeed clear in the 
Euclidean picture: for instance, particles coming from the past are rotated to come from the 
future, and can be re-interpreted as going towards the future if their charge is conjugated (a 
particle moving forward in time is the same as its anti-particle moving backward in time). 
In general, the parity invariance property (79) can be generalised to mixed states using the 
trace definition and (79) on the summand in the trace. 

Finally, we may also obtain the transformation property of twist fields. Since a a is 
spinless, and since <t q +i iQ and cr aiQ , +1 have opposite spins, it is clear that a a {x, y)<J a +\ (0) 
and <J a +i.a{x,y)o a ,a+i{fy transform under rotation by rotating the coordinates, without 
phase factors associated to the spin. There is an extra subtlety, however, having to do with 
the fact that we only ask for invariance of the state under rotation by the unique angle ■k. As 
a consequence, the correlation functions may depend on the direction of the cut emanating 
from the twist field, which we have chosen to be towards the right according to (21). After 
a 7r rotation, the cut goes in the opposite direction. This can be represented by the presence 
of an extra U(l) transformation operator, so that we have 

Ra a (x, y)tf = a a (-x, -y)e 2 * laQ (80) 

where Q is the Hermitian U(l) charge. The operator e 2ma< ^ is effectively a branch cut along 
the whole x-axis, cancelling out the twist fields cut to the right and reinstating the cut to 
the left. From (70) (and remembering that relation (74) inverts the positions of operators), 
this leads to 

Ra a+x , a {x,y)R) = e-^ a+1 ^a a+ha (-x,-y)e 2 ^ (81) 
Ra a>a+1 {x,y)R) = e^-Wacc+ii-x^y)****. (82) 
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TS = exp 



in 

y 



d9 (£>! [6)D\ {9) + D+ (0)£>_ {9)) 



Hence, a parity transformation gives us, using translation invariance of the state and the 
last relation (73), 

H(x,y) = (a a+ i tOl (x,y)a atOl+1 (0)) 

where the superscript c means that the state is as in the right-hand side of (79) (or its 
generalisation to mixed states). Likewise, 

G(x,y) = (a a (x,y)<r a+1 (0)) 

= (a a+1 (Q)<j a (-x,-y)e 4 ™ aQ ) c 

= (a a+1 {x iy )a a {Q)e^ a Qy. (83) 

Hence, parity invariance for these twist field correlation functions requires the following 
condition on the state: 

(... e AviaQ ) c =(■■■}. (84) 
If this condition holds, then wc find: 

G(-x,-y) = G(x,y), H(-x,-y) = H(x,y). (85) 

An example of a situation where there is translation and parity invariance is the finite- 
temperature state V = —H/T, with additionally the condition that 2a be an integer. A 
correlation function in a finitc-tcmperature state has the geometric meaning of a path integral 
evaluated on an infinitely long cylinder. The cuts of the twist fields extend towards the ends 
of the cylinder, and as such, affect the state representing the asymptotic values of the fields 
there (these asymptotic values are asymptotic states in the quantisation on the circle) . This 
means that correlation functions with cuts extending towards the right and the left are in 
general different. Hence, indeed parity symmetry imposes that there be no cut extending 
towards the right or the left. This imposes that the locality factor of the combined fields 
of the correlation function be one, hence that 2a be an integer. Note that by the relation 
between the Dirac theory and the Ising model, this gives the spin and disorder fields and 
their descendants. More generally, we see that a correlation function of the form (caUp), in 
a finite-temperature state and with the condition a + /3 G Z, is both translation and parity 
invariant. 

5.3 Ward identities 

Before embarking on our analysis of the Ward identities we first note that any correlation 
function of the form (cr a>Q! ±i (x, y)crp,p (0)) , (a a>a +i(x, 2/)cr a , a +i (0)) or (ct q , q _i (x, y)a a>a - 1 {Q)) 
vanishes, as a consequence of U(l) charge conservation. 

In order to derive a set of differential equations for the functions F(x, y), F(x, y), G(x, y) 



and H(x,y) we start with the following Ward identities, 

([Z,a%(x,y)aZ +1 (x,y)a* +l!a (0)al a+1 (0)})=0 (86a) 

([Z, [P,^(x,y)at +1 (x,yM +1<ol (0)at, a+1 (0)]) = (86b) 

([Z,[P,a^(x,y)at +1 (x,y)}a^ +lia (0)a% a+1 (0)])=0 (86c) 

([Z, [P, [P,**(x,y)a* +1 (x,y)]]al a+1 (0)a% +lta (0)]) = (86d) 

([ZAP,^(x,y)at +1 (x,y)][P,al a+1 (0)at +1 , a (0)]])=0 (86c) 
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and expand them using the action of Z obtained in Section 4.2. 

To ease our notation it will be useful to introduce the differential operator 

£>(/, 9) ■= \ {.fddg + gddf - dfdg - dgdf) (87) 

and also write T>(f) := T>(f, /). We note that 

V(e f ) = e 2f 8df. (88) 

From the Ward identities (86) we obtain, respectively, the equations 

e 2yTla H 2 + G 2 - FF = (89a) 

FdF - (dF)F = (89b) 

(dF)F-FdF = (89c) 

e 2 ™" (V(H) - 2m 2 H 2 ) - V(G) + V{F, F) = (89d) 

#***v{H)-V{G)-V{F,F) = 0. (89e) 

These look similar in structure to equations derived in [9] for the Ising spin-spin correlator 
except that here we have 4 functions to solve for. 

As F and F are only functions of x and y, Equations (89b) and (89c) can be solved to 
give 

F = k 2 F (90) 

where k is the constant of integration. This relation is somewhat surprising as it tells us the 
two correlation functions behave in the same way, at both short and long distances, despite 
the fields having different dimensions. Using the form factors for the twist fields k can be 
calculated exactly Inserting (18) between the fields in the correlation functions allows us to 
compare terms involving equal numbers of particles. Clearly the zero particle terms differ 
by the constant k 2 where 

k = C ^±±m 2a+1 (91) 
c a 

and it is a simple exercise to show that the same is true for all subsequent terms. We do not 
have yet a convincing intuitive explanation of the relation (90); it simply follows from our 
analytic-continuation definition of descendent twist fields. 

5.4 Differential equations for correlation functions 

Our next step towards an equation for the correlation functions is to define the functions % 
and (p via 

kF + G = e x cosh(<^) 

kF-G = e x sinh(<^) (92) 

so that (89a)becomes 

e 27ria H 2 = ie 2x sinh(2( /3 ). (93) 

It is now possible to write (89d) and (89e) in terms of \ an( i V only: 

sinh(2y>) (ddx + dipdip) + cosh(2<,c>) (ddip — coth.(2(p)d<pdtp) — m 2 sinh(2</j) = 094a) 
sinh(2^) (ddx + Zdtpdtp - ddip) + cosh(2(^) (ddtp - 2 coth(2^)9^ - ddx) = 094b) 
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and we can use (94a) to eliminate \ from (94b), leaving 

ddip - (1 + coth.{2y))dipdip + m 2 sinh 2 (2v3)(l - coth(2<p)) = 0. (95) 

In order to retrieve (3) from (95) it is necessary to eliminate the dcfidcj) term. To do this 
we note that 

/"(¥>) 



/'(¥>) 

admits the solution 



= 1 + coth(2</>) (96) 



— ni 

ddtp= — sinh(2V). (98) 



/(v>)=«r + lo 8 ^ i _ g4y + i j S 2^. (97) 
It is then straightforward to show that ip satisfies 

2 

To obtain (4) we set 

e E ° = 2A:F(a;, y) = e x (cosh(</?) + smh(tp)) (99) 

so that 

ddT, a = dd X + ddip. (100) 
Using (94a) and (95) it is possible to show that 

9 

— 77 1 

ddT, a = — (1 - cosh(2V)) (101) 

and so we have derived the known formula for the correlation functions (up to a constant 
factor which can be absorbed into the definitions of either the twist fields or the co-ordinates 
z and z). 

Following our method we notice that it is now possible to parametrise G(x, y) = (a a (x, y)a a -\ 
in terms of a partial differential equation involving ip. Setting 

e s ° = 2G(x, y) = e x (cosh((^) - smh(^)) (102) 

gives 

a5E' a = - 2dtpdip(l + coth(2^)) 
2tanh 2 (^) 



cosh(2-0) - 1 



dipdif). (103) 



This is a new result which is easily found following our construction. We also note that 
H(x,y) can be written in terms of S Q and Ti' a , via the algebraic relation (89a), and so it is 
also parametrised by the partial differential equations given above. 

6 Summary and perspectives 

As stated in the introduction, we have derived Equations (3), (4), (6) and (8). One natural 
extension of this work is the derivation of a differential equation for 

(a a (x,y)tT P (p,0)) for a£p,p±l 
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as is presented (for < a, /? < 1) in [2]. This calculation, however, is more complicated than 
that of the case a = (3, as some of the simplifications of Section 5 can not be applied. Some 
work is in progress in this direction [22]. 

Further, the derivation presented here is very direct as it involves only the symmetries of 
the model together with the action of the conserved charge Z on products of twist fields. It 
is also very general as it essentially only requires the density matrix to be an exponential of a 
quadratic form. It seems very likely that it could be applied to many other free field models 
possessing non-trivial twist fields. In particular, it would be very interesting to study, with 
such methods, twist fields associated to nonabelian symmetries. 

Also, the techniques used to calculated form factors for the descendent twist fields could 
be applied in other free field models, or in other quantum states. In particular, this may 
lead to ways of calculating vacuum expectation values of primary twist fields more generally, 
in direct analogy to our calculation of c a . Besides being useful again in the study of twist 
fields associated to nonabelian symmetries, it could be applied for instance in conjunction 
with the notion of finite-temperature form factors developed in [5], leading to expectation 
values at finite temperatures. 

Finally, a natural development is to determine the initial conditions that fix the solution to 
the differential equations corresponding to the correlation functions of interest. In particular, 
in analogy to what was done in [7] following [9] for the Ising model, one should connect the 
finite-temperature form factors of £7(1) twist fields in the Dirac theory to the initial scattering 
data (or Jost function) associated to the sinh-Gordon equation. 
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A Massless limit: CFT description and OPEs 



A.l Bosonisation and leading OPEs 

If the mass is set to zero in the equations of motion (9), then the fields ^>n and ^> R become 
holomorphic (right-moving), and the fields and ^ L anti-holomorphic (left- moving). In 
this case, the non-zero two-point functions of Dirac fields become exactly 



(vac|T ^ R (x 1 ,y 1 ) 1 $ R (x 2 ,y2) |vac) m=0 
(vac|T ^l{x 1 ,y 1 )^ L (x 2l y2) |vac) m=0 



z\ - z 2 
1 

Z\ - z 2 



(104) 



Multi-point functions are simply evaluated using Wick theorem. In particular, the right- and 
left-moving parts factorise. 

As is well known, it is fruitful to reproduce these correlation functions using bosonisation. 
We introduce holomorphic and anti-holomorphic operators e a {x,y) and e a (x,y), a £ K. 
They have bosonic statistics and the following correlation functions: 

(vac|e ai (a;i,2/i)---e Qn (a;„,2/„)|vac) m= o = ^o,£ jQj \\{-i{ z j ~ z k)T 3ak , 

j<k 

(vac|e ai (a;i,j/i)---e Q „(a:„,y„)|vac) m=0 = S Q ^, aj (i(zj - z k )) a]0lk (105) 



24 



for yi > y 2 > ■ ■ ■ > Dn, where the power functions are on their principal branches. These 
correlation functions imply the following operator product expansions: 

T[e a {x,y)e a ,{x',y')\ ~ (-i(z - z')) aa ' (l + -^{z - z')&\ e a+a ,(x' ,y') 

\ a + a' J 

T[e a (x,y)e a ,(x',y')} ~ (i(z - z')) aa ' ( 1 + -^—,{z - z')B') e a+a ,(x' ,y'), (106) 

\ a + a' J 

as well as the following equal-time exchange relations for the operators e a (x) and e a (x): 

e ai (xi)e a2 (x 2 ) = e^ a ^ s ^-^e a2 (x 2 )e ai ( Xl ), 

e ai (xi)e a2 (x 2 ) = e" a ^- s ^ x ^e a2 (x 2 )e ai ( Xl ) (107) 

in both cases for x\ =/= x 2 . 

The basic principle of bosonisation is the equivalence of such correlation functions with 
those of the Dirac fields for appropriate choices of as. More precisely, with the identifications 
V R = e-^/V^e-i, ^ R = e~ m / i Lo R e 1 , * L = e^Vc^i and = e^u^e-i, where oj r 
and u>l are pure phases, the correlation functions involving ^r, ^ r are reproduced, and those 
involving *l as well. In particular, the exchange relations (107) reproduce the fermionic 
statistics for exchanges in the right- and left-moving sector separately. However, in order to 
reproduce the correlation functions involving a mixture of right- and left-moving Dirac fields, 
we need to introduce extra elements guaranteeing fermionic exchange relations amongst right- 
and left-movers: so-called Klein factors. For this purpose, we will use quaternionic elements 
i , j and k , with both >f r and ^ R proportional to j , and both $l and W L proportional to i . 
In any given sectors, such quaternionic elements can be seen as being part of the pure phases 
wr,l, hence do not affect the correlation functions (using , etc.). For correlation 

functions of mixed sectors, the anti-commutation i j = — j i guarantees the correct statistics 
of Dirac fields. Since only even numbers of right- and left-moving Dirac fields arc involved 
in non-zero correlation functions, the quaternionic clement k never remains. The precise 
identification can then be written as follows: 

* fl = -e~ <,r/4 w~ 1 Je_i, V R . = e-^^ujeu $ L = e^wiei, tf* = -e"/ 4 ^ 1 ie_ x . 

(108) 

Here, we have used the U(l) invariance of the Dirac theory in order to remain with only one 
arbitrary phase uj. 

The bosonisation relations allow us to put in the same framework the [/(l)-twist fields 
discussed in Section 2. The basis for the identification is the comparison between the ex- 
change relations (21) and (107). Denoting by Q the [/(l)-charge of the Dirac theory, with 
[Q, ^ R L ] = — ty R L , [Q, $r,l] = ^r,l, we can simply make the identification: 

a a =e a e a e-" aQ . (109) 

The action of Q on the operators e a and e Q is given by 

e tcQ e a e~ lcQ = e- lca e a: e lcQ e a e~ lcQ = e tca e a . (110) 

Note that relation (109) along with the correlation functions (105) are in agreement with the 
normalisation (25). 

Correlation functions (105) can be realised using the Heisenberg algebra: one usually 
write e a (x,y) = e %aLf,R ^ x ^ and e a (x,y) = e~ tcapL ( x,v ^ where <pr,l are (appropriately nor- 
malised) holomorphic and anti-holomorphic free massless bosonic field, respectively. The 
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U(l) symmetry in terms of such fields corresponds to a shift of both tpn and ipi by the same 
constant, and the arbitrariness of the phase u in (108) can be thought of as an arbitrariness 
under shifts of tpn and ipi by opposite constants. Also, in terms of the bosonic fields, the 
[/(l)-twist fields is simply a a = e lQ v e - I7ra Q ; with tp := ip R — ip L . 

From relation (109) along with the correlation functions (105), it is possible to determine 
the leading coefficients of operator product expansions (OPEs) involving Dirac fields and the 
£7(1) twist field. Taking into consideration the exchange relations (21) and the correlation 
functions (105), we find 



T 



^ R (x,y)a a (0)] ~ Hz) Q (7 Q+lia (0), 



T[* R (x,y)<T a (p)] ~ (-iz)- a a a - lta {0), (111) 
where we define, for a € R, 

& a +i,a ■— e ui)e a+ ie a e ~', 

_ — i7r/4 —1- - —i-KaQ /-i -i r,\ 

<r a -i,a ■= -e ' oj }e a -ie a e w . (112) 

The power functions on the right-hand sides in (111) are on their principal branches; note 
that there is agreement between the phases occurring from exchange relations (21) and the 
phase differences occurring through the cuts of the principal branches. Further, we also find 



T 



^ L (x,y)a a (0)] ~ -iuj- 2 ie- i ^(iz)- a a a , a _ 1 (0) 



T[#L(x,y)a a (0)] ~ -M 2 ke™ Q {iz) a a a , a+1 (0). (113) 

Note that relation (33) can readily be obtained from (108), (112) and (109). However, in 
order to simplify the OPEs (113), we need to specify both k and uj. 

The CFT has both the U (1) symmetry, and a symmetry under constant shifts of <p. Hence, 
there does not seem to be any fundamental way of determining the phase uj in the CFT 
context. There are non-zero correlation functions involving odd numbers of both right- and 
left-movers, if twist fields are also inserted. Such correlation functions will involve remaining 
k factors, as well as remaining phases u>. Although fixing these correlation functions would 
partly lift these ambiguities, there is no clear principle as to how to fix them. 

In the massive theory, there is an (almost) unambiguous way of writing the element k 
explicitly as an operator. Indeed, k is an operator that anti-commutes with all Dirac fields, 
that commutes with the twist field a a , and that squares to —1, so we may simply write it as 

i=eie mQ (114) 

for some sign e = ±. In particular, we have (vac|k|vac) = ei. It is important that we can 
fix k , because there are unambiguous non-zero correlation functions involving odd numbers 
of both (massively deformed) right- and left-movers; in the bosonised language, these in- 
volve a remaining operator k , which needs to be well-defined. We note that in the present 
calculations, we do not need to determine the sign e. 

In addition to being able to fix k , in the massive theory, the symmetry under constant 
shifts of ip is broken. This symmetry breaking should allow us to evaluate the phase uj. The 
symmetry breaking is simple to see by writing explicitly the mass term in the Hamiltonian: 

— I dx (#Jj(a;)*£(x) - ^l(x)^ R (x)j = — — / dx cos{ip{x) + 2v) 
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with u =: e lv . In order to determine u, we simply define tp by specifying that the twist field 
a a has real and positive vacuum expectation value. Then, we use the fact that the field e %v 
has real and positive vacuum expectation value if and only if the integrand in the mass term 
has minima at tp = mod 2tt (so that the main contributions, in a path integral perspective, 
are at values of tp symmetrically distributed around 0: values of e lv add up to a positive 
real number). Hence, the integrand in the mass term must be costp and the coefficient in 
front of the integral must be negative in bosonic sectors (i.e. replacing k by (vac|k |vac)). 
Therefore, the requirement is that 

u? = -e. (115) 
Combining (114), (115) with (113), we find, independently of e, the OPEs 



T[*l(x,y)a a (0) 
T[* L (x,y)a a {0)] 



-(*2)- a <7 a , a _i(0) 
-(iz) a a a . a+1 {0). 



(116) 



It is a simple matter to use the definitions and relations found above in order to define 
other non-derivative descendent fields of the Dirac theory, and to establish the leading and 
first sub-leading terms in the OPEs with the Fermion fields. We define recursively <T a ± n ,a 
for integer n > via 



T 



{-iz) 



' a+n-j- 1,q 



(0), 



This gives the expressions, for aeR and n > integer, 

_ / — ITT / 4 ?\n - — iirotQ 

cr a +n,a ■= (e ' Wj) e a+n e a e v , 



(117) 



(118) 



(these are in agreement with (109) and (112) in the cases n = and n = 1 respectively). 
From these expressions along with the rules (108), (114), (115) and (110), and with the OPEs 
(106), we find the OPEs (for a - /3 G Z): 



T 



T[* R (x,y)<T aiP (0)] 



H) e(«-/»-l/2) ( _ ij8) -« ^ + ^zd) a a . u 



(0) 



T 



^ L (x,y)a ail3 (0) 



T[^ L (x,y)(r a ,p(0)] 



( _ i) e(a-^-l/2) ( _ 1) a-^+ 1( ^^ 



1-/3 



zd a a ^-i(0) 



1 + /3 



f5U ai/}+ i(0XH9) 



where 6 is the Heaviside step-function, 8(7) = 1 (7 > 0), (7 < 0), and the derivatives on 
the right-hand side are with respect to the arguments of the twist fields, which are afterwards 
put to 0. In these OPEs, the next correction term is equivalent to a term of order z 2 (for 
the first two) or z 2 (for the last two) inside the large parenthesis (but generally involve 
non-derivative descendents of the twist fields written) . 

We may write these OPEs more generally by putting the twist fields at position x', y' and 
by replacing z by z — z' and z by z — z! in the power functions on the right-hand side. From 
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this, taking derivatives with respect to z' and z', we obtain OPE's with first-derivatives of 
twist fields: 



T 



0)0^(0)1 = {-i) s v- a - 1 '%a)(-iz) a - 1 *^ 1 , p (o) + o(ir +2 ) 



T y)fcr aij9 (0)] = (-i) e C- a, - 1 / 2 5(-l) a -^ +1 (i J 8)(«)- /, -V a ^_i(0) + 0(z^+ 2 ) 

T [* L (x,y)d* a AQ)] = (-O e(a "^ 1/2) (-i) a "^ +1 (-^)(^)^V a ,0 +1 (o) + 0(*"+ 2 ). 

(120) 

Note that the first non-trivial corrections to the leading OPEs, corresponding to pure deriva- 
tive of twist fields, are exactly zero. 

A. 2 Action of the charge Z 

The OPEs above along with the integral expression (59) of the charge Z allow us to evaluate 
the action of Z on the product of twist fields belonging to different copies as follows. First, 
in the CFT the charge Z decomposes into a sum of two conserved charges, Z = Zr + Zl- 
where Zr^l involves only right-moving and left-moving Fcrmions respectively. Second, by 
the usual principles of QFT, commutators of local charges with local fields can be evaluated 
through contour integrals and the residue theorem. Using the replacements — iz > s and 
iz n- s respectively, the commutators with Zr and Zl are recast into the following contour 
integrals: 

[Z R ,O(0)] ^ ^Ids(*R(x,y)$ R (x,y) + y R (x,y)$ R (x,y))o(0) 

[Z L ,O(0)] ^ -^Ids(^L(x,y)^ L (x,y) + ^ L (x,y)^L(x,y))o(0). 

The OPEs (119) above immediately give rise to Equations (62), and the OPEs (120) provide 
actions of Z on products of twist fields involving one derivative that agree with those written 
in Section 4.2. 

B Functional equations without translation and parity 
symmetry 

In this appendix we write down the differential equations arising from the Ward Identities 
(86) in the absence of translation and parity symmetry The loss of these symmetries means 
that we no longer have all the nice symmetry properties discussed in section 5.2. As a result 
we must introduce a more appropriate, but cumbersome, notation. We define 

F a S p(.xi,yi,X2,y2) ■= (<r a ,0(xi,yi)<T y ,s(x2,y2)) (121) 

and first note that the lack of translation symmetry prevents us from always placing the 
second field at (0, 0), and so the correlation functions must remain functions of four variables. 
It should also be noted that 

F t,a{ x UVl,%2,V2) = (<7a,a(xi,yi)<rp,p(x2,y2)) = (<Tp,p(X2,y2)<rcc,a(xi,yi) 

) ^F^(x 1 ,y 1 ,x 2 ,y 2 ) (122) 
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and so the Ward identities can not be expected to simplify as much as they did in the 
previous case. 

The Ward identities in this case are 



([ Z , O-a.aOl, 2/l )^If l,a (^2 , 2/2^* Q + 1 (^2 , 2/2)]) 


= 


(123a) 


([Z, [^^* Q (^i ; yi)ffa+i, Q +i(^i^i)] CT * + i,Q(^2,2/2)CT* Q+ i(a;2,2/2)]) 


= 


(123b) 


([Z, [^^* Q (xi,yi)ff* +l Q+1 (zi,yi)]CT* +1 Q (x 2 ,2/2)CT* Q+ i(a;2,2/2)]) 


= 


(123c) 


^ [-Pi [-P, a O&l, 2/1)^+1^+1 (^1 , 2/l)]]0-* o+ i (»2 , 2/2)cr*+i : « (x 2 ,2/2)]) 


= 


(123d) 


^ 0-*a(^l, S/O^+l.o+l (^15 yi)][-P\ ^a+l^, 2/2)0-a+l,a(2 ; 2, 2/2)]]) 


= 


(123e) 



and their expansions in terms of correlation functions are cumbersome. To ease our notation 
the dependence on (x\, yi, x<i, y%) will be implicit in all F^'t, as operators at (xi,yi) arc 
always to the left of those at (22,2/2)- We also introduce the symmetric, bilinear differential 
operator 

T> it j(u, v) = udidjV + vdidjU — diudjV — djudiV (124) 
where Zj = -?r{xj + iyj) and 

a i = — and di = —. (125) 

Now we can write (123) as 

pa+l,a cia,a+l _ pa+l,a+l na,a . pa,a pa+l,a+l n fine; \ 

r o,«+l r a+l,a Oi,ct r otLotl t o.a o+la+1 ~ U \, lzUc V 



a p«+l,o pa,a+l 
^Zl^ajtt+l r a + l,a 



— F, 



a+l,a 
a,a+ 



1 "2i r ci+l,i< 



a pia+l.aH 



1 p a , a 
r a+l,cH 



™ + l,a+la 



!+l,Q!+l 



pa, a a pO+l,a 
^ct.a L 'zi- r a+l,a- 



fl = (126b) 



} pa+l,a pO,a+l 
'z\ r a ,a+l r a+l,a 



F 



a+l,a 
a,a+ 



1 "2i r a+l,i< 



TP 
l r < 



+ F a ' Q 

o pa+l,a+l i?ot,a 
y zi r a+l,Q + l r a,a 



a+l,aH 
a, a 



F a,a 



+ l t '2i r Ct, 



-l.Q+l 



a+l,a+l 



1^,^=0 (126c) 



7-) /p«>" + l T7ia+l,a\ 2 j-ia,a+l p a+l,a 

^l.l^a+l.a ) f a ,a+l ) ~ m *a+l,a ^a.a+1 

©i,i(*S£. F:+l;:tl) + ^(FS% ha+ \^i, a+ i) = (126d) 

T>i,2(F£#,F2$%) ~ V lt2 (F%Z, F^tD ~ ^{FS% ha+ \F^, a+x ) = 0. (126e) 

It should be noted that this system is under defined as there are 6 independent functions, 
each of two variables, and only 5 equations. It is possible to find more Ward identities leading 
to non-trivial relations between the correlation functions but as we are unable to identify any 
integrable structure in these equations we only include the above equations for completeness. 
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